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ABSTRACT 

The  known  differential  equation,  which  describes  the  spreading  of  a  circular 
electron  beam  in  an  accelerating  Held  under  space-charge  repulsion,  is  de¬ 
rived  in  a  simplified  form  so  that  the  general  solution  can  be  precomputed. 
The  solution  is  discussed  and  is  given  is  a  diagram.  All  auxiliary  calcula¬ 
tions  are  listed  in  order  that  the  beam  radius  can  be  found  for  any  set  of 
initial  conditions.  Only  the  diagrams  and  a  slide  rule  are  needed.  Another 
independent  method  of  computation  for  the  same  purpose,  which  may  be  use¬ 
ful  under  special  conditions,  is  also  included. 
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COMPUTATION  OF  THE  SPREADING  OF  AN  ELECTRON  BEAM 
UNDER  ACCELERATION  AND  SPACE-CHARGE  REPULSION 


I.  INTRODUCTION 

In  electron  optical  devices  in  which  there  is  a  high  electron  current  density,  as  in  cathode-ray 
tubes,  the  field  due  to  the  electrons  themselves,  i.e.,  space  charge,  is  an  important  factor  in 
the  determination  of  the  shape  of  any  electron  beam.  The  diverging  effect  of  space  charge  in  a 
region  which  is  field  free  is  a  standard  presentation  in  many  books  on  electron  optics.  The  equa¬ 
tion  of  motion  is  usually  attributed  to  Watson  (cf.  Hef.  3) 

The  space  charge  effect  in  an  accelerated  electron  beam  is  also  important  but  does  not  seem 
to  have  been  treated  extensively  in  the  literature.  We  will  investigate  a  circular  electron  beam 
in  an  accelerating  field.  The  electrons  have  a  certain  initial  velocity,  and  enter  the  axial  elec¬ 
trical  field  -  for  example,  through  a  hole  in  a  plane  electrode.  The  electrons  are  now  subjected 
to  constant  acceleration  and  mutual  repulsion.  The  equation  of  motion  has  been  derived  and  solved 
by  H.  Moss'  who  gives  several  examples  for  the  conditions  found  in  cathode-ray  tubes 

Such  an  example  is  defined  by  a  set  of  six  independent  parameters  By  Moss'  metnod,  every 
example  is  solved  by  computing  a  function  numerically  and  then  integrating  it.  For  a  given  beam 
length,  the  limits  of  integration  i  .ust  be  found  by  trial  and  error  The  mathematical  difficulties 
have  apparently  preven'ed  any  wider  application  of  this  solution. 

In  the  course  of  voltage  breakdown  studies  we  became  interested  in  field  emission  as  an  ini¬ 
tiating  mechanism  The  high  current  density  of  field  emission  makes  the  space-charge  effect 
especially  important.  In  this  connection  it  was  desirable  to  compute  many  cases  which  would 
have  been  laborious  by  Moss'  method  A  simplified  method  has  been  found  for  computing  the 
shape  of  an  accelerated  circular  electron  beam  with  space  charge,  and  we  believe  that  it  may  be 
useful  for  other  applications  This  method  is  presented  here  in  detail. 

In  Sec  11  the  mathematical  derivation  and  a  discussion  of  the  resolving  function  are  given. 
Section  III  covers  the  numerical  computation  of  the  resolving  function.  The  function  is  repre¬ 
sented  in  Figs  2  and  3, 

The  reader  who  wants  to  solve  a  specific  problem  may  omit  Secs.  II  and  Ill  and  simply  follow 
the  instructions  given  in  Sec.  IV  where  all  necessary  steps  for  the  application  of  the  resolving 
function  are  listed.  A  numerical  example  is  also  given  in  Sec.  IV. 

In  Sec.  V  a  step-by-step  procedure  is  described  which  is  based  upon  the  formula  for  a  beam 
in  a  field-free  space.  This  is  a  second  independent  solution  of  the  same  problem. 
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n.  MATHEMATICAL  ANALYSIS 
A.  Equation  oi  Motion 

For  the  sake  of  completeness,  we  will  derive  our  simplified  procedure  from  the  fundamentals, 
following  closely  the  original  deductions  of  H.  Moss.  After  our  Eq.  (12)  we  deviate  from  Moss' 
analysis 


Fie-  I  •  Gaomstry  for  ipocs-chorga  spreading  of  a  boom 
in  an  •Uclric  field. 


We  look  at  a  thin  slice  (thickness  Az)  of  a  circular  electron  beam  (Fig  1)  and  assume  that 
in  any  plane  perpendicular  to  the  beam  axis  the  electrons  arc  monoenergctic.  The  total  electro¬ 
static  flux  due  to  the  space  charge  in  the  slice  is 


y 


I)  (is 


o 


lAz 

7. 


(1) 


where 


I)  flux  (displacement  density), 

S  surface  of  0  cylindrical  slice, 

i>  charge  within  the  slice, 

I  beam  current, 
z  velocity  of  the  beam  electrons. 

Since  we  assume  that  the  field  strength  and  the  flux  are  constant  in  axial  direction,  the  end  discs 
do  not  contribute  and  F.q  (1)  becomes 

D  2irrAz  -  ^  (2) 


The  radial  electrical  field  is  then 


F  -J2.  =  _L  I 

r>rad  e  c  TrrrJ. 
o  o 


and  the  radial  force  on  the  electrons  is 

el 


F  =  eE 


rad  e  2*rz 


(3) 


(4) 
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The  equation!)  of  motion  of  the  electrons  are  with  an  external  axial  field  E 


d2r 

dt2 


d2z 


dt 


We  integrate  Eq.  (6) 


t„2»rz  m 
o 


i  =  —  Et  +  z„ 
m  o 


(5) 


(6) 


(7) 


z  =  -  |  t2  t  i  t  =  A  V  *  Bt 

m2  o  2 


Equation  (7)  inserted  in  Eq.  (5)  yields 

el 


d2r  .  _ _ 

dt2  eo2»rm  Et  +  zo) 


(8) 


(9) 


This  is  Eq.  (6)  of  H.  Moss1  except  for  a  factor  1/4  »cu  due  to  the  use  of  different  units.  In  order 
to  simplify  the  differential  equation,  we  introduce  the  abbreviations 
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K  = 


el 


e  2»m 
o 


A  =  ^  E 
m 


131 

1*1 


(10) 


(11) 


B  =  io  ■ 

and  the  dimensionless  variables 
A 


I— I 

1  sec' 


R  =  r 


T  =  t 


'/BK 

A_ 

B  • 


Z  =  z 


F 


We  insert  Eqs.  (13),  (14)  and  (14a)  into  Eqs.  (8)  and  (9)  and  obtain,  respectively. 


Z  =  |  T2  +  T 


(12) 

(13) 

(14) 

(14a) 

(Mb) 
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and 


d^R  1 

t2  '  T  +  1 


(15) 


The  dimensionless  Eq.  (15)  could  be  solved  as  it  stands  by  some  numerical  method,  but  it  would 
still  give  a  two  parametric  set  of  curves,  which  is  difficult  to  compute  and  to  represent. 

B.  Formal  Solution 

We  write  instead 


(T  +  1) 


3/2  d2R 
dT2 


=(  — M'1 

'  Vt  +  if 


(16) 


The  solution  of  this  Eq.  (16)  is  given  in  the  literature,  e  g.,  by  Kamke.2  We  may  call  this  solu¬ 
tion  a  formal  one  since  it  is  given  in  terms  that  have  no  direct  physical  meaning  We  introduce 
the  new  variable  into  Eq.  (15) 


W  =  _R_  =  r  ,  (17) 

n/tT1  *  B 

and  obtain  from  Eq.  (IS)  after  one  integration 

(T  +  I)2  W2  =  |  W2  +  2  J  dW  =  j  W2  «•  2  In  W  t  C  (18) 


Separation  of  variables  yields 

P  *  dW 

\  * - ™ - -  ln(T  ♦  1)  ♦  C.  .  (19) 

j  w  n  z  ■ 

0  w  +  2  lnW  +  c 

The  integral  represents  our  resolving  function.  Before  we  compute  it  numerically,  we  will 
make  it  independent  of  the  parameter  WQ  [Eq.  (28)|.  The  boundary  conditions  are 

t  =  0  ,  r=r  .  r  =  r_  =  r '  i  , 

o  o  o  o  ’ 

T  =  0  ,  R  =  R  ,  R  =  R 

o'  o 

A  dot  on  a  capital  letter  means  differentiation  by  T;  a  dot  on  a  lower  case  letter  means 
differentiation  by  t. 

The  boundary  conditions  inserted  in  Eqs.  (18)  and  (19)  give 


C  =  W  2  -  ~  W  2  -  2  In  W 

O  4  O  O 

(20) 

o 

it 

O 

(20a) 
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where 


W  =  R  -  i-  H  r  n/B/K  -  4-  r 
o  a  i  o  o  '  To  sjjfc 


(21) 


W  -  K  =  r  --- 

°  °  °  VBK 


(21a) 


C.  Discussion  of  the  Resolving  Function 

The  different  possible  signs  in  Eq.  (19)  require  special  care.  We  rewrite  Eq.  (18) 

*  W2  +  2  in  W  t  C 

W  '  — ^ — f— j -  .  (22) 

The  denominator  (T  +  1)  is  always  positive.  WQ  is  known  from  the  boundary  conditions  and  de¬ 
termines  which  sign  must  he  applied  at  the  start  of  the  computation  WQ  also  indicates  the  direc¬ 
tion  in  which  we  have  to  integrate  since  T  in  Eq.  (19)  must  be  positive  and  can  only  increase. 
Therefore, 

if  W  <  0,  the  negative  sign  in  Eq.  (19)  applies  and  W  <  WQ, 
if  Wo  >  0,  the  positive  sign  in  Eq  (» 9)  applies  and  W  >  WQ. 

It  remains  to  be  investigated  whether  at  any  later  time  T  there  may  be  a  change  of  the  sign  of 
W,  i.e.,  the  point  W  =  0  needs  special  attention.  We  introduce  the  index  n  and  write 

Wn  0  (23) 


and  also  hold  the  equations  ready  for  use 
R  1  R 


W 


vT  +  1 
Ft 


W 


(T  +  l)3^ 
R„ 


"  J T  +  1  2  (T  +  1)3/2 

V  n  n 


=  0 


We  insert  Eq.  (23)  into  Eq.  (22)  and  conclude  that  for  Wn  =  0  either 


(24) 


(25) 


T  = 


or 


^  W  2  +  2  In  Wn  +  C  =  0  (26) 

The  last  case  is  the  interesting  one;  it  means  that  at  Wn  the  integrand  of  Eq.  (19)  goes  to  infinity. 
An  investigation  of  the  limit  of  the  integral  in  Eq.  (19)  is  given  in  Eq.  (31),  and  it  is  shown  that 
the  integral  plotted  vs  W  remains  finite  and  has  only  an  infinite  slope  at  W  =  Wn. 

At  this  point,  W  -  Wn,  we  have  to  reconsider  the  question  of  signs  in  Eq.  (19).  It  is  physi¬ 
cally  evident  that  R  is  always  positive  and  increases  monotonically  or  first  decreases  to  a 
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minimum  and  then  increases  '  With  this  in  mind,  it  can  be  seen  from  Eq.  (24)  that  if  W  started 
negative  and  became  zero  at  Wn>  it  can  only  be  positive  afterward.  That  means  that  from  WR 
on,  the  positive  sign  applies  in  Eq.  (19),  and  the  direction  of  integration  is  reversed  so  that  W 
increases  from  this  point  on. 

With  W  determined  by  Eq.  (26)  we  can  define  the  integral  independent  of  the  parameter  WQ, 
with  only  the  parameter  C  left.  We  write 

f(W)  =  -  ■  1 -  ,  <27> 

W2  +  2  In  W  +  C 


c w  t 

F(W)  =  \  . riT—r 

Wn  W2  4  2  lnW  + 


dW 


In  order  to  compute  F(W)  from  f(W)  it  is  necessary  to  discuss  (he  limit 


lim  f(W) 
W-W 


(28) 


because  at  this  point  the  integrand  becomes  infinite 

t 


f(W) 


s/t  (Wn 


♦  tf  ♦  2  ln(Wn  ♦  c)  *  C 


Vf  <Wn  +  Tc2+2,n,‘ 

since  1/4  WZ  +  l  lnW  +  C  -  0  according  to  Eq.  (26), 


lim  f(W)  = 
e-0 


(l/4)-(l/Wn)Z 
(Wn/2)  +  (2/Wn) 


1  - 


il'W 


±  1 

4-  W2 


-2-  +  — 
2  '  W, 


n  '  J 


(29) 


*  No  cronovsr  )•  possible  under  the  assumed  conditions. 
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The  second  term  in  the  expression  is  kept  only  to  estimate  the  error.  We  insert  Eq.  (29)  into 
Eq.  (28) 


F(Wn  ♦  .)  f(Wn 


+  c)  d« 


— D - 


We  take  as  an  approximation 


F(Wn  +  c )  =  2  V7 


for  which  the  relative  error  will  be  smaller  than 


(30) 


(31) 


A  = 


(32) 


If  a  maximum  error  A  is  allowed  (and  known),  Eq.  (32)  gives  the  e  up  to  which  we  apply 
Eq.  (31),  i.e.,  Eq.  (28)  becomes 


F(W)  =  2  <€ 


i 


r 

jw  _• 


dW 


JR 


+  2  InW  +  C 


(33) 


This  equation  is  the  basis  for  the  numerical  computation  of  F(W)  as  it  is  explained  below. 

It  is  worthwhile  to  note  that  the  minimum  value  of  W,  which  we  have  called  W_,  does  not 

n 

correspond  to  the  minimum  of  the  beam  radius  r  or  R.  The  minimum  beam  radius  must  sat¬ 
isfy  the  equation 


OR  -  R  -  0 
dT  -  R  -  0 


With  Eq.  (24)  we  can  write 
R 


W 


1  R 

2  (T  +  l)3/2 


W<Rmin>  =  -i 


min 


<Tmin  +  1> 


372 


7 


and  with  Eqs,  (17)  and  (22), 


T  +  1 


W  <  W„  , 
o 

which  determines  W(Rmin). 

ID.  NUMERICAL  COMPUTATION  OF  THE  RESOLVING  FUNCTION 

We  have  chosen  a  set  of  values  for  the  parameter  C  between  the  limits  -20  and  +20.  This 

covers  all  practical  applications  we  can  anticipate.  For  every  C  we  first  compute  W  by  solving 

"  „  3 

Eq.  (26)  with  the  Ncwton-Kaphson  method.  The  residual  error  is  made  smaller  than  10  per¬ 
cent.  Next  we  determine  the  initial  step  width  c  from  Eq  (32)  assuming  A  -  10~*  and  check 

whether  c  <  W  which  is  the  condition  for  Eq.  (32)  being  valid.  If  t  >  W  we  take  r*  =  0.9  W 
n  "  n  n 

instead  as  initial  step  width.  Now  we  compute  F(W)  with  a  recursion  formula: 
from  Eq.  (31 ) 


W,R  )-  *W(,W 

",Kmln  2  T  +  1 


This  reduces  to  the  equation 


21nW<Rmin>  +  C  =  0 


Ft  2 


yW  , 


(34a) 


from  Eq.  (29) 


t. 


•J7  [W 


/w - - 

-H  +  -2. 

7  2  Wn 


(34b) 


f,  = 


(35  a) 


+  2  InWj  +  C 


wi  =  wi-l +  * 

> 

(35b) 

Fi  ■ 

Fi-1 + 1 

tfi-l + 1  ffi 

if  i  is  even 

(36) 

Fi-1  -  I 

fi-2  +  l  ffi-l  +  3  <fi 

if  i  is  odd 

(37) 

Equation  (36)  is  a  simple  extension  of  the  integral  over  the  next  step  with  the  trapezoidal  rule; 
Eq.  (37)  has  the  following  appearance  if  we  insert  successively  the  preceding  equations  for 
i-1,  i  -  2,  etc. 


Fl  +  !«fl*T 


+  f  ef,  ♦...+ 


i'h. 


i + j  <fi 


(38) 
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This  is  Simpson's  Rule.  After  every  second  step  (i  odd)  we  may  adjust  the  step  width  t.  Then 
we  have  instead  of  Eq.  (38) 


Ff  =  Fj  +  j  tjfj  +  ]  €  1^2  +  3  {i^3  +  3  £3^3  +  3  €  3^4  +  3  f3^5 


+  +  5  £i-4fi-2  +  T  £i-2fi-2  +  I  £i-2fi-l  +  I  ei-2fi 


(39) 


This  is  again  Simpson's  Rule  with  varying  step  width. 

Equations  (34)  through  (37)  have  been  coded  for  the  IBM  7090.  The  computed  F(W)  is  plotted 
in  Figs.  2  and  3. 

Once  F(W,  C)  is  known  and  a  set  of  initial  conditions  is  given,  we  apply  the  formal  solution 
as  follows.  We  insert  Eqs.  (20a)  and  (28)  into  Eq.  (19) 

F(W.C)-  F(Wo,C)  =  ln(T  +  1)  ,  (40) 


and  invert  Eq.  (14b)  to 

T  +  1  =  +  sTTTTZ  (41) 

We  insert  Eq.  (41)  into  Eq.  (40)  and  write 

F(W.C)  =  lnt/TTIZ  +  F(Wo>C)  (42) 


Since  W(),  C  and  Z  are  defined  by  the  initial  conditions  and  F  is  known,  Eq.  (42)  yields  W.  We 
write  Eq.  (17)  with  Eq.  (41)  inserted  as 

K  =  W  t/TTl  =  W  <A  +  2Z  (43) 


In  this  equation  we  insert  W  and  find  R.  By  inserting  Eqs.  (13)  and  (21a)  into  Eq.  (43)  we  obtain 
finally 


4  2Z 


(44) 


We  can  now  demonstrate  the  step-by-step  application  of  the  formal  solution  starting  with  specific 
initial  conditions 

IV.  APPLICATION  OF  THE  FORMAL  SOLUTION  TO  A  PHYSICAL  PROBLEM 

To  start  the  computation,  we  must  know  all  the  parameters  which  define  the  beam  at  the 
starting  point.  These  are: 

Beam  current 
Initial  beam  radius 

Initial  beam  divergence 

Initial  potential  of  electrons 
Accelerating  electric  field 


dr, 


U 

E. 


9 


to 


v»  W  for  C  nigotivi. 
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In  order  to  find  the  beam  radius  r  at  a  given  beam  length  z  we  compute  successively  the 
following  quantities: 


i  =  = 

o  v  m 


A  e  /£n2,rm 

VV  -  r  —  -  r  —  F  /  —2 - 

°  °  “m  ‘N  *oel 


|Eq.  (21a)] 


K *  {'0  •/! =  ro 


•  j  , 

z  el 


C  K  -  W  It  -  2  In  W 

O  O  O  ( 


(Fig  4)  , 


A  K  — 

Z  ;  A-  ,  z  _JL  "1 

*  u2  Z  .2 
B 


IEq.  (14a)] 


|F(Wo.C)|  is  taken  from  Fig.  2  or  Fig.  3.  The  sign  is  still  indeterminate; 
we  must  take  the  sign  of  W  . 


W  -  It  -  t  W 
o  o  2  o 


[Eq.  (Zl)] 


W 


F(Wo.C)  -  ^  ■  |F(Wo.C)| 


(45) 


F(W,  C)  F(Wo,  C)  +  In  •JTTJZ  [Fig.  5  and  Eq.  (42)] 


W(F)  is  taken  from  Fig.  2  or  Fig.  3 

iL  =  W  y,  +  2Z  [Eq.  (44)  ] 

o  o 

bet  us  look  at  two  examples.  In  the  first,  the  initial  condition  correspond  to  the  focusing  of  a 
post -acceleration  CRT;  in  the  second,  to  field  emission. 
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Initial  Conditions 


Example  1 

Example  2 

I  =  25  •  10‘6amp 

1  ■  10'6  amp 

r„  =  1  ■  to-3 m 

1.785  ■  10’9i 

ro  r~m 

0 

U  =  1  •  103 V 

5  v 

Ez  =  5  •  103  v/m 

5  •  107  v/m 

z  -  0.25m 

e* 

O 

1 

3 

Parameters 

Example  l7 

z  =  1.873  •  107  m/sec 
o 

W  ^  0.7223 
o 

R  =  -  1.435 
o 

C  =  3.781 
Z  =  0.628 

|F(Wo,C)|  =  0.5100 

W  -  -1.796 
o 

F{W  ,  C)  =  -0.5100 
o 

lrWFVTZ  =  0.4080 
F(W,  C)  =  -0.1030 
W(F)  =  0.1791 

Results 

Example  1  Example  2 

--  0.3040  845 

The  result  In  Example  1  is  compared  in  Fig.  6  with  that  of  H.  Moss  for  the  same  conditions. 


Example  Z 
1.236  •  10^  m/sec 

2.43  •  10*1 
0 

2.83 

5  •  103 
0 

0 

4.62 

4.62 

20.5 


t  These  numerical  values  have  mere  significant  digits  than  one  could  read  from  the  diagrams,  and  the  result  Is 
more  precise  then  usually  needed. 
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Ffg.  6.  Relative  beam  radius  vt  length 
for  condition*  of  Example  I.  the  curve 
It  that  of  H.  Meu.  The  point  Indicator 
our  remit  for  Z  *>  230  mm. 


V.  ANOTHER  INDEPENDENT  WAY  TO  COMPUTE  BEAM  SPREADING 


Instead  of  supposing  a  homogeneous  field,  we  now  assume  many  slices  perpendicular  to  the 
axis  of  the  beam,  where  the  potential  is  constant  within  every  slice.  At  the  cut  the  potential 
jumps  to  the  potential  of  the  next  slice.  Within  each  slice  the  beam  spreads  according  to  the 
formula  for  a  field-free  space.  At  every  cut  we  have  to  apply  the  law  of  electron  optical  refrac¬ 
tion  because  the  angle  of  divergence  is  reduced  by  the  jump  of  the  potential. 

The  appropriate  formula  for  beam  spreading  is  taken  from  a  recent  r&sumfe  of  Ivey  5  and 
after  correction  of  an  obvious  misprint  it  reads 


1.021 


(46) 


where 


r  =  minimum  beam  radius  at  the  "waist", 
m 

S  =  axial  distance  from  the  "waist”  of  the  beam, 
r  =  beam  radius, 

V  -  potential  in  kv, 

1  «  beam  current  in  amperes. 


.  _  .  ("1  lnfr/r  )  2 

<p(V)  *  l(for  our  purpose). 


The  function  F(r/rm)  is  taken  from  Jahnke-Emde4  and  plotted  in  Fig.  7.  We  represent  it  within 
the  range  r/rm  =  1. .  .  2  by  the  simpler  expression 


(47) 


which  is  also  plotted  in  Fig.  7. 

With  F*  inserted,  Eq.  (46)  becomes 


2.14 


1/2 


(4S) 


or 


0.220 


S‘  I 

71  V172 

m 


+  1 


The  divergence  of  the  beam  is 


dr 

da 


0.44 


I 

^72 


(49) 


(50) 


Since  for  our  conditions  of  field  emission  the  spreading  of  the  beam  depends  only  upon  the  current 
density,  we  may  choose  a  very  small  initial  diameter  and  a  correspondingly  small  current;  then 
we  have  a  very  small  angle  of  divergence  For  a  sufficiently  small  beam  divergence,  we  can 
compute  the  radius  at  the  end  of  the  it4>  slice  by  adding  separate  components: 


(a)  The  radius  at  the  entrance  of  the  slice, 

(b)  The  increment  due  to  space-charge  repulsion  within  the  slice,  computed 
for  an  initially  parallel  beam, 

(c)  The  increment  due  to  the  initial  beam  divergence  which  equals  the  beam 


Thus  we  obtain 


r.  --  0.22 


»i  =  0  44 


7^- 

(ASj)2 

ri-t 

I 

V.3/2 

_ 

4 

♦  r.  -  ♦ 


ri-t  V/‘ 


. dr .  .  „  (i7 1 

<d^>i-t  *Si  J  —  • 

dr,  /ST 

ds’i-i  J  Vj 


(51) 


(52) 
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Assuming  initially  (dr/dg)  *  0,  we  can  compute  r  step  by  step.  We  have  only  to  watch  that 

ri  <  2ri_1  because  otherwise  our  representation  of  F(r/rm)  fails.  It  is  best  to  start  with  a  very 

small  AS,  e.g.,  10'5  mm  and  to  double  the  stepwidth  every  cycle.  In  such  a  way  for  E  =  50kv/mm, 

z  =  1  mm.  V_  =  5ev,  r_  =  1.78S  •  to'*1  mm,  I  *  to*6  amp  we  obtain 
o  o 

This  happens  to  be  just  the  value  that  we  have  found  with  the  graphical  method  in  Example  2 
(Sec.  IV).  (In  general  we  have  to  allow  for  a  reading  error). 

The  coding  of  Eqs.  (51 )  and  (52)  for  the  IBM  7090  is  reproduced  in  the  appendix. 

VI.  CONCLUSION 

The  two  examples  shown,  and  many  more  which  have  been  computed,  indicate  that  both 
methods  give  the  same  results  within  the  limits  of  the  reading  error. 

The  first  method,  which  uses  the  resolving  function,  is  appropriate  if  one  wants  to  quickly 
compute  a  few  examples.  For  the  exploration  of  a  variety  of  conditions,  the  first  method  also 
offers  a  fairly  simple  way  to  obtain  one  result  at  a  time,  which  is  useful  for  the  choice  of  the 
next  input.  Of  course,  the  first  method  is  always  recommended  if  no  electronic  computer  is  at 
hand. 

The  second  method  automatically  gives  the  complete  beam  shape  up  to  the  length  under  con¬ 
sideration,  but  it  always  involves  the  use  of  a  high  speed  computer.  It  is  therefore  most  suitable 
for  longer  series  of  examples  as  they  occur  by  systematic  variation  of  parameters  in  small  steps. 
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APPENDIX 

FORTRAN  PROGRAM  FOR  SOLUTION  OF  EQS.  (51)  AND  (52). 


ELECTRON  BEAM 

C  DIMENSION  DELSOI2000).  SOtfOOO),  V(2000),  ROP(2000),  RP(200), 

CAPRPI2000),  CAPRI2000),  RHOtfOOO),  RHOP(2000).  FIRHO<2000), 
FLOGtfO),  FLSS(2000) 

16  READ  INPUT  TAPE  2.  1,  SO(l),  DELSO(i).  RHO(l),  FI.  U,  RHOP(l), 

SEND  WRITE  OUTPUT  TAPE  1,  3 
FM=0  434294 
PI-3.14159265 
1  =  1 

11  1=1+1 

DEL.SO(I)=2.0+  DEl.SO(l-l) 

12  SO0)=SO(l-l)+DELSO(I) 

V(I)-(SO!I)+SO(I-l))*  U/2.0 

ROP(I)=Fl/SQRTF(V(I>*  *  3)*  (DELSO(I)/RHO(l-l))*  *  2*  0.220+1.0 
IF(ROP(1)-2.0)  10,11.  11 

11  DELSO(I)=DEI.SO(1)/2.0 

GO  TO  12 

10  RP<I)=0  44*  DELSO(I)*  Fl/(RHO(l-l  )*  SQRTFIV(I)*  *  3)) 

C'APRP(I)=RHOP(I-l  )*  SQRTF(V(1-1  )/V(I)) 

C'APR(1)=CAPRP(1)*  DELSO(I) 

RHO(I)=ROP(l)*  RHOU-1  )+CAPR(I) 

RUOP(I)=RP|I)+CAPRP(I) 

FIRHO(I)=FI/(PI*  RHO(I)*  *  2) 

FI.OG(I)=LOGF(FIRHO(I»*  FM*  0.2 
FLSS(I)=LOGF(SO(I)-SO(l»*  FM*  0.2 

WRITE  OUTPUT  TAPE  3,  2.  SO(I),  RHO(I).  FIRHO(I),  FLOG(I),  FLSS(l) 
WRITE  OUTPUT  TAPE  14,  4,  Fl.SS(l),  FLOG(I) 

IFfI-2000)  14,  15, 15 

14  IF(SOd)-SEND)  13. 15,  15 

15  GO  TO  16 


V  FORMAT  STATEMENTS 

1  FORMAT  (5E10.4.2F6.2) 

2  FORMAT  (3E20.6.2F20.6) 

3  FORMAT  (10H  S16X,  4H  RH014X,  6H  l/RHOSX,  12H  LOGFU/RHO)  10X,  1 

10H  LOG(S-SO)/) 

4  FORMAT  (2F20.6) 

END 
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EXPLANATION  OK  THE  SYMBOLS  USED  IN  THE  FORTRAN  CODE 

Fortran  Equations  (51)  and  (52) 

SO(l)  SQ 

SEND  SK 


SO(l) 


s 


i-1 


DELSO(l) 

V(l) 

ROP(I) 

RP(I) 

CAPKP(I) 


RHO(I) 


r 


i-1 


RHOP(I) 


'ds'i-1 


FIRHO(I) 

FLAG(l) 

FLSS(l) 


I 


0.2  lotJ0  — 

*ri-l 

0.2  log  jo  (Sj.j-S^ 


for  plotting  purposes 


The  last  three  symbols  are  not  explained  in  Eqs.  (51)  and  (52).  They  are 
only  used  to  plot  the  results  in  a  convenient  form. 
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